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A Brownian particle moving across a porous membrane subject to an oscillating force exhibits 
stochastic resonance with properties which strongly depend on the geometry of the confining cavities 
on the two sides of the membrane. Such a manifestation of stochastic resonance requires neither 
energetic nor entropic barriers, and can thus be regarded as a purely geometric effect. The magnitude 
of this effect is sensitive to the geometry of both the cavities and the pores, thus leading to distinctive 
optimal synchronization conditions. 

PACS numbers: 05.40.-a, 05.10.Gg 



Stochastic resonance (SR) is by now a textbook ex- 
ample of how noise can best enhance the response of a 
bistable system to an external drive [T]. Historically the 
research on SR focused mostly on systems with purely 
energetic potentials, either continuous or discrete. How- 
ever, as pointed out in Ref. [2] , in soft condensed matter 
and in a variety of biological systems [3], particles are 
often confined to constrained geometries, such as inter- 
stices, pores, or channels, whose size and shape can affect 
the SR mechanism 0]. Indeed, smooth confining geome- 
tries can be modeled as entropic (i.e., noise or tempera- 
ture dependent) potentials [5], capable of influencing the 
response of the system to an external driving force (see, 
for a review, Ref. [5]). 

Let us consider a Brownian particle freely diffusing in a 
two (2D) or three dimensional suspension fluid contained 
in two symmetric cavities connected by a narrow pore; it 
can switch cavity only by overcoming the entropic barrier 
determined by the geometric constriction associated with 
the pore. This is true even for ideal reflecting boundaries, 
that is, in the absence of an intrinsic energetic barrier. 
Burada et al. [I] have shown that entropic barriers signif- 
icantly contribute to the magnitude of the SR effect that 
occurs when a periodic force drives the particle across the 
pore; hence the term of "entropic SR" coined in Ref. [4j. 
The evidence reported there, however, hints at an inter- 
play of entropic and energetic barriers, rather than to a 
mere entropic effect. Indeed, SR was demonstrated there 
only under the explicit condition that the applied field 
of force had a dc component at an angle with the pore 
axis, say orthogonal to it, so that the particle tended to 
sojourn preferably against one side of the cavities. In the 
absence of such an additional symmetry-breaking force, 
no SR was observed. 

This remark raises the issue whether a bistable effec- 
tive potential is a necessary condition for SR to occur 
altogether [7j. In this paper we show that a Brownian 
particle confined to two distinct cavities divided by a 
porous medium, say a membrane, does undergo SR when 
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FIG. 1: (Color online) (a) Brownian particle confined to a 2D 
box divided in two compartments by a partition with an open- 
ing at the center, (b) Three different pore geometries used in 
our simulations: funnel, hole, and spout (left to right). A and 
xo denote, respectively, the cross-section and the thickness of 
the pore. 



driven by an ac force perpendicular to the membrane, 
even in the absence of external gradients and/or interac- 
tions with the walls (besides bouncing from the walls). 
At variance with ordinary SR, optimal synchronization 
between drive and particle oscillations for an appropri- 
ate noise level, only occurs at driving frequencies (am- 
plitudes) lower (higher) than a certain onset threshold. 
Moreover, such a manifestation of SR in higher dimen- 
sions requires adopting extremely sharp geometrical con- 
strictions to separate the two cavities, something akin to 
the pores obtained by puncturing a thin membrane. The 
magnitude and conditions of the effect reported here are 
sensitive to both the geometry of the cavities and the 
cross-section of the pores, thus allowing a direct control 
of the synchronization mechanism. 

The overdamped dynamics of a Brownian particle in 
2D is modeled by the Langevin equation 



df/dt = -A(t) e x + VD £(t) 



(1) 



where e x , e y are the unit vectors along the x, y axes and 
= (^x{t),^y(t)) are zero- mean, white Gaussian noises 
with autocorrelation functions (£i(i)£j(f')) = 2<Sy<J(t — t') 
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FIG. 2: (Color online) Geometric stochastic resonance, (a) 
x(D) versus D for different values of the ac drive parameters 
A at fi = 0.01 in (a) and fi at A = 0.045 in (b). Other 
parameters are: s_l = j/z, = 1 and A = 0.1. The dashed 
curve in (a) represents our predicted asymptotic decay x/Aq 
for D — >• oo. In (b) we display our predictions for: the range 
of variability of x(0) (horizontal arrows), the SR peak posi- 
tion D max (vertical arrows), the SR peak height 5;(-D ma x) (top 
line), and the decay law x(D — > oo) (dashed line); see text. 



with i,j = x,y. Equation ([I]) has been numerically in- 
tegrated for the two-cavity container sketched in Fig. [T] 
with reflecting walls [8] and a single opening, with differ- 
ent geometries, placed at the center of the partition wall. 
In the presence of an ac drive, A(t) = y4 cos(£1£), the 
Brownian trajectories embed a persistent harmonic com- 
ponent, x(D) cos[f2£ — <fr(D)}, whose amplitude, x, and 
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phase, (j), are plotted versus D in Figs. 

For simplicity, we start modeling the pore as a struc- 
tureless hole pierced in a zero-thickness wall. The occur- 
rence of a SR phenomenon is clear, albeit with some dis- 
tinctive features, as seen in Figs.[2]and[3| (1) x{D) peaks 
for an appropriate noise intensity, Z? max , with -D max a 
manifestly increasing function of O; (2) SR is restricted 
to A) > A c , Fig.^a), and Q < ti c , Fig.^b). This is an 
important difference with respect to ordinary SR, where 
there exist no such onset thresholds in the drive param- 
eters space; (3) The curves x(D) decay like -D -1 , Fig. |3j 
that is faster than in any ID bistable potential and more 
in line with SR in a discrete two-state model [3]. 



All properties listed above can be explained by simple 
geometrical considerations. We start noticing that when 
the ac force A(t) presses the particle against the walls 
of the container opposite to the dividing wall, then the 
average particle displacement (x(t)) approaches a square 
waveform with amplitude xl- Correspondingly, the par- 
ticle gets pushed against the container partition twice per 
cycle; if it goes through the opening, (x(t)} traces a sym- 
metric, two-sided square wave; if it doesn't, its average 
displacement is restricted to an asymmetric, one-sided 
square wave, on either the positive or the negative side 
of the dividing wall. On taking the Fourier series of (x(t)) 
with period 2ir/n, the amplitude of its fundamental har- 
monic component turns out to be, respectively, Axl/tt 
for the two-sided waveform, and 2x^ / it for the one-sided 
waveform. From this remark there follow immediately an 
upper bound to x, 



x < (4/7r) x Ll 
and the onset condition for SR, 

A /n > (4/tt) x L , 



(2) 



(3) 



Ao/fi being the driven oscillation amplitude of an un- 
constrained Brownian particle. The SR peak, x(D max ), 
approaches the upper bound in Eq. |2 1 at vanishingly low 
fi, see Fig. [2jb); whereas from Eq. (3|, for our simula- 
tion parameters one obtains A c = 0.013 in Fig.|2ja) and 
il c = 0.035 in Fig. [2]^b). The consistence of these analyt- 
ical results with the simulations is quite satisfactory. 

At variance with SR in a bistable potential, the parti- 
cle oscillations are not drastically suppressed in the zero- 
noise limit. An estimate for x(0) = Huid^qX^D) can be 
obtained by noticing that the probability for the particle 
to cross the pore and, therefore, to execute a full oscilla- 
tion is A/ul, whereas the probability to get trapped on 
either side of the partition is (1 — A/j/l). Accordingly, 

x(0) = (x L /2)(l + A/y L ) K (n), (4) 

with k — 4:/tt for O — )■ , see Eq. |2|, and n = 1 for 
fl — > f2 c . This defines a relatively narrow variability 
range for 5;(0) as a function of ft. Well above the onset 
threshold namely for fl <C fi c , for D — > the curves 
in Fig. ^h) clearly tend to the upper bound of a;(0) in 
Eq. Q; correspondingly, the upper bound of noise am- 
plification via SR is x(D)/x(0) < 0.5. 

To check these predictions we modified the container 
geometry. In Fig.[3]the response amplitude x(D), in units 
of xl, is plotted versus D for different A's, panel (a), and 
Xl's, panel (b). In agreement with Eq. Q, x(0)/xl in 
panel (a) grows linearly with A, while in panel (b) is 
seemingly insensitive to xl- Note that x(0) approaches 
the upper bound ([2| only for A — > 1, thus implying that 
geometric SR happens for any finite pore width. 

The dependence of the decaying tails of x{D) on A 
and Xl is also consistent with our geometric interpre- 
tation. Indeed, the D — > oo behavior of x{D) can be 
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FIG. 3: (Color online) Geometry dependence of this SR. x(D) 
versus D for different values of A (a) and xl (b). Other 
parameters: t/l = 1, Ao = 0.045, Q = 0.005 and (a) xl = 1, 
(b) A = 0. 1. In both panels 5;(0) is very close to our prediction 
of Eq. Q with = 4/tt. 

analyzed in terms of the average time, t±(D), for an 
unbiased Brownian particle to diffuse across one com- 
partment. Such time constant is easily obtainable by 
analytical means [8], T\(D) — x 2 L /iD. In the presence 
of strong noise, the geometric constriction exerted by 
the pore grows ineffective; diffusion along the x-axis is 
then described by an ac forced damped Brownian mo- 
tion with effective damping constant r-j" , so that the 
corresponding VL component of (x(t)) gets suppressed 
both in amplitude, A Tx(D)/y/l + [SIti(D)] 2 , and phase, 
4>{D) = arctan[Ori(-D)]. For Q <^fl c , this leads to the Q- 
independent estimate, x(D) ~ AqT\{D), plotted in both 
panels of Fig. [2] (dashed curves). Note that n is inde- 
pendent of A and quadratic in xl, so that for D — > oo 
the tails of the curves of Fig. [3]Ja) collapse onto one curve 
and those of Fig.[3^b) scale proportional to xl- 

The SR peaks of x(D) occur at a certain value of the 
noise strength, Z3 max , see Fig. [2] which weakly depends 
on Aq and grows with fl. This property distinguishes ge- 
ometric SR from ordinary SR, as here the x(t) transitions 
between stable states (cavities) are not regulated by an 
activation rate over an energetic barrier of the Arrhenius 
type. In the absence of a drive, Aq — 0, the relaxation is 
characterized by some exit time, t(A), which is necessar- 
ily inverse proportional to D. Optimal synchronization 
between pore crossings and external drive requires that 
the particle switches compartment twice during one A(t) 
cycle pQ, namely 

r(A) = tt/O. (5) 

As a consequence, -D max is proportional to f2. 

For a quantitative analysis, we propose to define t(A) 
as the mean time a Brownian particle, uniformly dis- 
tributed in one compartment, first crosses the pore. This 
quantity has been numerically computed and plotted in 
Fig. |4|a) for Aq = and different container geometries. 
On inserting the corresponding numerical data for r(A) 
into Eq. (J5|, we predicted the values of D max marked in 
Fig. [2^b) by vertical arrows. 



The dependence of r(A) on the parameters A and Xl 
is also instructive. For A = we recovered the ID 
limit, n, as it should. At very small pore cross-sections, 
t(A) turned out to scale like (y^/A) 1 / 2 , no matter what 
the compartment aspect ratio xl/ul- However, the ratio 
xl/ul controls the interplay between diffusion along the 
x and y axis. For long compartments, xl/ul 3> 1, the 
exit time is of the order of t\ until quite low A, whereas 
for narrow compartments, xl/i/l <C 1, the exit time is 
dominated by the vertical diffusion of the particle to- 
wards the mid-section of the cavity, where the pore is 
located, i.e., t(A) ~ (j/l — A) 2 /12D for A not too close 
to and 1 [also shown in Fig. Qa)]. To this regard, we 
remark that such a dependence of the exit time ouxl/ul 
is not inconsistent per se with a reduced ID description of 
the Brownian diffusion along the x-axis. However, none 
of the assumptions introduced in Ref. [3] to accommo- 
date for SR in the framework of the Fick-Jacobs kinetics 
|10j . apply to the present case. In particular, no en- 
tropic barrier can be naturally attributed to the pore. 
The mechanism investigated here is rather an irreducible 
2D geometric effect, for which the term geometric SR is 
more appropriate. 

The phase delay, 4>{D), of the fundamental harmonic 
component of (x{t)) with respect to A(t) is plotted in the 
insets of Fig. [2j By comparing the corresponding curves 
for x{D) and 4>(D), we notice that: (1) <fi(D) also ex- 
hibit SR-like peaks, which shift to higher D on increas- 
ing f2 [Fig. [2jb)] , while being almost insensitive to Aq 
[Fig. [2ja)]. Note that in a two-state model 4>(D) would 
be a monotonically decreasing function of D [1] . (2) The 
phase peaks get suppressed at low Q, where the SR sig- 
nature of x(D) is the most prominent; (3) In the zero 
noise limit, contrary to x(0), <fi(0) strongly depends on fl 
and Aq; for D — > oo the expected asymptotic behavior, 
4>{D) ~ firi(-D), was recovered (not shown). The non- 
monotonic behavior of <fr(D) signals the appearance of 
two-sided oscillating trajectories in the averaging ensem- 
ble of (x(t)), with the particle crossing more and more 
frequently the pore in unison with A(t). According to 
our interpretation of the phenomenon under study, phase 
peaks are natural SR precursors. However, if the SR on- 
set condition (|3| is not met, 4>{D) tends to 7r/2, as ex- 
pected, as the particle is restricted to oscillate inside one 
cavity most of the time pQ. 

In real experiments at small-length scales, the geome- 
try of the partition wall and its opening(s) are often not 
fully controllable. In Fig. |4]jb) we plot x(D) versus D for 
the three different pore shapes sketched in Fig. [TJb). In 
order to enhance geometry effects, we simulated a square 
waveform A(t) also with amplitude Aq and angular fre- 
quency Q. Moreover, the pore width, xq, was taken not 
too small, lest the different pore geometries become indis- 
tinguishable. Funnel-like pores tend to channel the Brow- 
nian trajectories through the opening, no matter what 
their cross-section; thus, no SR evidence was detected, 
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FIG. 4: (Color online) Dependence of SR on the pore ge- 
ometry, (a) t(A) in units of ri for different cavity lengths. 
The drive has been switched off and the initial position of 
the particle uniformly distributed within one cavity; averages 
were taken over 2 x 10 trajectories. The dashed curve repre- 
sents the vertical diffusion time, (i/l — A) 2 /12D (see text), for 
the shortest cavity. Inset: scaling law ti/t(A) oc (j/l/A) 1 ^ 2 
for A — > 0. (b) x{D) versus D for different pore geome- 
tries: spout (circles), hole (squares), and funnel (triangles). 
A(t) is a square waveform with Q = 0.01 and amplitude as 
in the legend. Inset: x(D) versus D for a spout-like pore at 
Ao — 0.045 and different fi. Other simulation parameters are: 
x jj — y L — 1, A = 0.1, and x = 0.1 (where it applies). 

as anticipated in Ref. [4]. On the contrary, spout-like 
pores exhibit enhanced SR peaks and a distinct depen- 
dence on the cavity geometry. Indeed, under the pressure 
exerted by the drive, the exit time through a spout of 
width xq is of Arrhenius type, namely, a function of the 
pore cross-section, proportional to e A o x °/ D . Owing to 
the SR condition |5]), Anax is now expected to grow al- 
most linearly with A and logarithmically decreases with 
fi, in qualitative agreement with our simulations. Note 
that, as soon as the pore shape comes into play, i.e., for 
xqAq 3> D, the Aq dependence of the exit time cannot 
be ignored any more and a quantitative analysis of these 
results requires going beyond the approximations of the 
linear response theory [T]. In the opposite limit, the pore 
can be well modeled by a simple hole, as we initially did. 

We expect that geometric SR can be best demon- 
strated in vortex superconducting devices [11] . This class 
of artificial devices is presently attracting growing in- 



terest because of potential applications to flux qubits, 
SQUIDs and superconducting rf filters. Ion-beam tech- 
niques allow to fabricate a superconducting sample with 
two vortex boxes connected by a thin pore of almost any 
geometry. Vortices are trapped inside the boxes with 
binding energy of the order of ^Lt/X 2 , where $o is 
the magnetic flux quantum, A is the London penetra- 
tion depth, and L t is the depth of the two vortex traps. 
The vortex density n = _ff/$ j is controlled by the inten- 
sity H of the applied magnetic field. In the dilute limit, 
H < <&o/A 2 , the vortex- vortex interactions become negli- 
gible, so that the transport properties of a single trapped 
vortex are not overshadowed by many-body effects, ac 
drives and noise sources can be easily implemented as 
Lorentz forces generated by independent electric currents 
injected into the sample parallel and perpendicular to the 
pore axis. Detection of SR under such experimental con- 
ditions is regulated by the applied current sources only; 
in particular, the noise parameter D can be varied inde- 
pendently of the constant operating sample temperature. 
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